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Abstract. We consider the generation of large-scale magnetic fields in slow-roll 
inflation. The inflaton field is described in a supergravity framework where the 
conformal invariance of the electromagnetic field is generically and naturally broken. 
For each class of inflationary scenarios, we determine the functional dependence of the 
gauge coupling that is consistent with the observations on the magnetic field strength 
at various astrophysical scales and, at the same time, avoid a back-reaction problem. 
Then, we study whether the required coupling functions can naturally emerge in well- 
motivated, possibly string-inspired, models. We argue that this is non trivial and can 
be realized only for a restricted class of scenarios. This includes power-law inflation 
where the inflaton field is interpreted as a modulus. However, this scenario seems to be 
consistent only if the energy scale of inflation is low and the reheating stage prolonged. 
Another reasonable possibility appears to be small field models since no back-reaction 
problem is present in this case but, unfortunately, the corresponding scenario cannot 
be justified in a stringy framework. Finally, large field models do not lead to sensible 
model building. 



PACS numbers: 98.80.Cq, 98.70.Vc 
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1. Introduction 

Magnetic fields are present on various scales in the Universe from planets, stars, galaxies, 
to clusters of galaxies, see e.g., Refs. [U (2J [3] for reviews. In galaxies of all types, 
magnetic fields with the field strength ~ 10~ 6 G, ordered on 1 — 10 kpc scale, have been 
detected [3, H] . There is also some evidence that they exist in galaxies at cosmo logical 
distances [5]. Furthermore, in recent years, magnetic fields in clusters of galaxies have 
been observed by means of the Faraday rotation measurements (RMs) of polarized 
electromagnetic radiation passing through an ionized medium [HI [TJ [8] . Unfortunately, 
however, RMs give only the product of the field strength along the line of sight and the 
coherence scale, and we cannot find them independently. In this situation, the strength 
and the scale are estimated on 1CT 7 — 1CT 6 G and 10 kpc— 1 Mpc, respectively. It is 
interesting that magnetic fields in clusters of galaxies are as strong as galactic ones and 
that the coherence scale may be as large as ~ Mpc. 

Since the conductivity of the Universe through most of its history is large, the 
magnetic field B evolves conserving magnetic flux as B oc a~ 2 , where a(t) is the scale 
factor. On the other hand, the average cosmic energy density p evolves as p oc a~ 3 . 
Hence, the scaling of the magnetic field is given by B oc p 2 / 3 and this relation should 
also be understood as giving the value of B at different spatial locations characterized 
by different energy densities. The present ratio of interstellar medium density in 
galaxies p ga i to p and that of inter-cluster medium density in clusters of galaxies p c \ are 
Pgai/p — 10 5 — 10 6 and p c \/p — 10 2 — 10 3 , respectively. Consequently, from these relations, 
we see that, using the law mentioned before, namely B = -B ga i,ci(Pgai,ci/p)~ 2//3 , the 
required strength of the cosmic magnetic field at the structure formation, adiabatically 
rescaled to the present time, is 10~ 10 — 10~ 9 G in order to explain the observed fields 
mentioned before. On the other hand, in general, if the galactic dynamo mechanism is at 
play, then seed fields with a present strength of only 10 -22 — 10 _16 G is required [9l fT0lfTT] . 

The question addressed in this article is that of the origin of these magnetic fields 
and whether there exists simple, physically well-motivated mechanisms, which could 
provide the target values of B just described above. Generation mechanisms of the seed 
magnetic fields fall into two broad categories. One is astrophysical processes and the 
other is cosmological physical processes in the early Universe. The former exploits the 
difference in mobility between electrons and ions. This difference can lead to electric 
currents and hence magnetic fields. If the scale of cluster magnetic fields is as large 
as ~ Mpc, however, it is unlikely that the origin of such a large-scale magnetic field 
is in astrophysical processes and we are tempted to assign the origin of cosmological 
magnetic field to the processes during inflation in the early Universe which stretches 
coherent scales exponentially [121 H31 HH HSl HSl HZl HH [HI [20]. 

During inflation quantum scalar (density) perturbations [ZD, E2J [231 El] and tensor 
(gravitational waves) perturbations [251 126] are generated and stretched to cosmological 
scales by subsequent accelerated expansion. Usually, however, large scale fluctuations 
in electromagnetic fields are not generated, because they are conformally invariant [27J . 
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Hence, in order to generate cosmological magnetic fields during inflation, one must break 
the conformal invariance of the theory in an appropriate manner. A number of models 
have been proposed in this context so far [28l [29l [30l EH E21 [33l EH [35l [36l EZl EHl EHl SO] . 

Here we propose a simple model making use of the arbitrariness of the gauge kinetic 
function in supergravity [41] . Indeed, the general supergravity Lagrangian including 
vector super-multiplets contains the following term 



where W a a is the gauge field strength superfield with spinor index a and gauge group 
index a. As already mentioned, the quantity f a b{&) is an arbitrary function of the chiral 
superfield $ [12] ■ This implies that the conformal invariance is naturally and generically 
broken as is most conveniently shown when the above term is expressed in terms of the 
superfields components 



Since we are interested in the production of magnetic fields, we consider a U(l) gauge 
field and, therefore, the gauge field strength do no longer carry a gauge index. As a 
consequence, there is also only one function /(</>) left. Moreover, in order to consider 
the minimalist case with the least possible assumptions, we naturally identify the scalar 
part (ft of the chiral superfield $ to the inflaton field. 

A related study was carried out in Refs. [29J, ED] but only in the particular case where 
f(4>) oc e^. Here, we study the most general case (at least, in the framework of single 
field inflation) and determine which form f(4>) should take in order to obtain a power 
spectrum of the magnetic field in agreement with the currently available observational 
constraints on large scales. Then, for each category of single field inflationary scenarios, 
we explore whether there is a natural model in supergravity (possibly string-inspired) 
that could reproduce the required shape of the gauge coupling and, at the same time, 
leads to the correct form of the inflaton potential. 

This paper is organized as follows. In Sec. [21 we describe our model action and 
derive the equations of motion from it. Then, the gauge field is quantized in the 
Coulomb gauge and the vacuum energy density of the magnetic field at the end of 
inflation is determined. The evolution of the magnetic energy density through the 
subsequent stages of evolution, preheating, radiation and matter dominated eras is also 
computed. In a next step, for each single field slow-roll inflationary models, using a 
simple ansatz depending only on one free parameter, we calculate the required shape 
of the gauge coupling function /(</>). In Sec. [31 we use the currently available limits on 
the magnetic field to constraint the free parameter mentioned above, that is to say to 
constraint the shape of the gauge coupling function. We also investigate the consistency 
of the model and determine under which conditions this one does not suffer from a 
back-reaction problem. In Sec. HJ we try to embed in particle physics the models that 
have been shown to be phenomenologically viable. Finally, in Sec. [5j we present our 
conclusions. Throughout this article, we use units in which c = h = Mq — ks — 1. 




(1) 




(2) 
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We also put the reduced Planck mass, Mq = 1/ V8ttG, to unity when it simplifies the 
notation. In this case, the full Planck scale takes m P1 = G -1 / 2 = y/8n. 

2. Reverse-Engineering of the Magnetic Field 



2.1. Magnetic Field during Inflation 

During inflation, following the considerations presented in Sec. [TJ the action of the 
system, namely a scalar inflaton field plus a U(l) gauge field, can be written as 



s [0, 4J 



d 4 X A 



gfsTP (0) F m F v 



■13 > 



(3) 



with F^ a = d^A a — d a A^, A a being the vector potential whose dimension is that of 
a mass [A^\ = M. The inflaton potential V((f>) and the dimensionless gauge coupling 
/(</>) are, a priori, arbitrary functions. Notice that we have changed the definition of 
the gauge coupling (now f 2 appears in the action rather than /) for future convenience. 
The equations of motion following from the above action can be written as 



(4) 
(5) 



% [V=gg' u g aP f (0) f vP ] = o . 

In the following, we consider that the electromagnetic field is a perturbatively small 
quantity (a "test" field) compared with the inflaton counterpart. This means that one 
can neglect the right hand side of the Klein-Gordon equation and work in a spatially 
flat Friedmann-Lemaitre-Robertson- Walker (FLRW) space-time, 

ds 2 = g^dx^dx u = -dt 2 + a 2 {t)dx 2 = a 2 {n) (-dry 2 + da; 2 ) , (6) 

where t is the cosmic time and rj the conformal one, since the gauge field being negligible 
there is no preferred direction. Therefore, the equations of motion are just the standard 
ones, namely 



H 



\tf + v(4>) 



a \ 1 

(/ / 3 

with Ti = a' /a and, for the inflaton field 

dV 



- 2 n 2 



+ 3H(t> + 



d<p 



dV 

, 0" + 2H<p' + — = . 



(7) 



(8) 



In order to simplify the calculations, we adopt the Coulomb gauge where A (t,x) = 
and djA>(t,x) = 0. In this case, one obtains 



Ai(t, x) + [ H + 2j ] A(t, x) - dj&Aitt, x) = 



A': + 2 f -A\ 



a 2 did 3 A, 







(9) 
(10) 
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where, of course, <9 J = g^ k dk = /a 2 )dk- If one defines Aj = fAj, then one can absorb 
the part proportional to the first derivative of the vector potential and put the equation 
of motion under the standard form of an equation for an oscillator. One gets 

A?-^-A i -a 2 d j d i A i = 0. (11) 

We see that, if / = 1, the equation of motion of A 7 reduces to that of an harmonic 
oscillator. Due to the conformal coupling of the gauge sector, the dynamics is non 
trivial only if / is not a constant. 

Let us now (covariantly) define the electric and magnetic fields seen by an observer 
characterized by the four-velocity vector u M . One has [13] 

Efj, = u v F^ v , Bp = -e liuK F VK , (12) 

where the tensor e Mi , K is defined by the relation 

£fu/K = Vhvk\ u ) (13) 
V^vkX being the totally antisymmetric permutation tensor of space-time with rj 0123 = 
(— g)~ 1 ^ 2 or r/0123 = ( — g) 1 ^ 2 - Therefore, for a comoving observer with = (1,0) (in 
cosmic time), one gets 

Ei = —Ai , Bi = -tijkdjA k , (14) 

(X 

with ei23 = 1. A remark is in order concerning the formula giving the magnetic field. 
Contrary to the previous expressions, it is not written in a covariant way and, therefore, 
the position of the indices is irrelevant even if implicit summation is still present. So, for 
instance, this formula simply means B\ = (^2^3 — <9 3 v4 2 )/a. As expected, the dimension 
of the magnetic field is two, [Bi] = M 2 . 

We now turn to the quantization of the system. It follows from Eq. (J3j that the 
momentum conjugate to the gauge field Ai(t, x) is given by 

7T%X) = ^- = /W(*)<MM) = f^i^A'^x) , (15) 

5Ai 

which also implies that iTi(t,x) = f 2 ((p)a 3 (t)Ai(t,x). For quantization, we impose the 
canonical commutation relation between A % (t,x) and 7Tj(t,x): 

[A% x), 7Tj(i, y)]=ij ^re*^ (V, - 6 jt ^j (16) 

= iSf^ (x-y), (17) 

where try^ is the transverse delta function introduced in order to have a consistent 
quantization in the Coulomb gauge. 

In order to define the creation and annihilation operators, one has to Fourier expand 
the vector potential operator. For this purpose, one introduces an orthonormal basis in 
space-time defined by (in conformal time) 

4=(? )- £ -(4)> ^Kt)- (18) 
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where A = 1,2 and where, by definition, o~ij£\£ J x = 1 (no summation on A). In the 
above equations, k l = ak l phys is the comoving wavenumber satisfying b id k % y = k 2 . In 
addition, e x e^ — e\ki = which expresses the Coulomb gauge in Fourier space. Notice 
that, without the factor 1/a in the previous definitions, the above vectors would not be 
correctly normalized in curved space-time. Then, one has the covariant completeness 
relation J2lZo dss^s 6 ^ = 9^ where g ss is just a factor which stands for +1 or — 1. 
Projected on space, the completeness relation reduces to 

j^e\{k)e 3X {k)+5^ = V 3 . (19) 



k 2 

A=l 

Finally, the expression of Ai(t, x) can be written as 

+ b[(k)A*( V ,k)e-* k - x ], (20) 

in terms of the annihilation and creation operators, b\(k) and b x (k), with k being the 
comoving wavenumber. The quantity e^(fe) is the transverse polarization vector and 
has been defined before and the time-dependent Fourier amplitude obeys the following 
equation of motion 



A"( V ,k)+(k 2 -ipj A(r),k) = 



(21) 



where A(f],k) = a(rj)f(jj)A(rj,k). The additional factor a in the previous definition 
originates from the presence of the vector €i\(k) in the Fourier expansion (!20l of the 
gauge field. The dimension of the Fourier amplitude of the gauge field is [A] = M~ x l 2 
since [AJ = M and [&a(&)] — M 3 ^ 2 (because the delta function, see below, has dimension 
equal to the inverse of its argument). In order to satisfy Eq. (|16|) . the creation and 
annihilation operators must satisfy 

b x {k),b{,(k')] =6 3 (k-k')5 xx ,, [b x {k),b y {k')} = [b{{k),b{,{k')\ =0, 

(22) 

where we have used Eq. (1191) . The normalization condition for the time-dependent 
amplitude A(t, k) should be chosen such that 

A(t, k)A*(t, k) - A(t, k)A*(t, k) = -i- . (23) 

f 2 a A 

Let us now calculate the energy density of the magnetic field. The stress energy 
tensor is given by 

T» = - -f(*)**W» ~ ^W-V*^ , (24) 

and the energy density associated to the "magnetic" part can be expressed as (of course, 
there are also contributions depending on d Ai, that is to say depending on the electric 
field, that we do not consider here) 

Too = W" (djA e - dtAi) (diA k - d k AA = f -^a 2 B,B» . (25) 
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As expected, one has [Tqq] = M 4 thanks to [Ai] = M. Using the Fourier expansion of 
the potential vector, straightforward manipulations lead to 

(T B °o) = ~ P f { l y J ^k\A{nM 2 ^9 mi 

2 

x ^2 ( eexk i ~~ e i xk ^ ^ mxki ~ e ^ k m) ■ (26) 

A=l 

The final step is to notice that 

2 

g i Jg mi ^ (e^kj - e jX k e ) (e mX ki - e iX k m ) 

A=l 

=4t^-^ k 't^<-4' (27) 

A=l A=l 

where we have used the space component of the completeness relation. The energy 
density, p B = — (T BO o), is therefore given by 

A{r],k) " 



i r+°° au 



2vr 2 J k 
1 f + °° &k 



a(rj) 

or, considering only the energy density stored at a given scale, 

lfM'>' k) = ^JM lAMl2 - (29) 

As expected, the fourth power of the physical wave number, /c p h ys = k/a(j]) appears in 
the above expression. 

Let us now compute the power spectrum in the case where the scale factor is given 
by a power law of the conformal time, 

1+/3 



a(rj) = a Q 



(30) 



rj_ 

\Vo 

The case (3 = —2 corresponds to de Sitter space-time. One can also show [44] that 
the spectral index of density perturbations is given by n s = 2f3 + 5. If we adopt the 
very conservative bound \n s — 1| < 0.1 then this amounts to choosing the index (3 such 
that —2.05 < j3 < —2. Notice that the above assumption for the scale factor is not 
as restrictive as it may seem at first sight. Indeed, the above law is also valid in the 
slow-roll approximation where one has [H] 

a(r]) oc l^r 1 " 61 , (31) 

where e\ is the first slow-roll parameter and is constant at first order. Therefore, our 
ansatz allows us to treat the case of slow-roll inflation as, for instance, large and small 
field inflation or even hybrid inflation in the inflationary valley. 



Generation of Large-Scale Magnetic Fields in Single-Field Inflation 



8 



Let us now discuss the form of the gauge coupling function /(0) that will be 
considered in this article. One assumes that 

f(v) oc a a , (32) 

where a is a free index. The motivation for such a choice is as follows. Unfortunately, the 
precise form of the magnetic spectrum is not known experimentally. Only upper limits on 
the amplitude at given scales have been obtained so far, see below. Therefore, one cannot 
take the function p B and reverse-engineer it exactly in order to find the corresponding 
f{4>)- On the other hand, it will be demonstrated below that the coupling function with 
a = 2 leads to a scale- invariant power spectrum, see also Ref. [39J. Clearly, since the 
amplitude of the primordial magnetic field is not strongly peaked over a specific range 
of scales (either small or large scales), a scale-invariant spectrum satisfies the currently 
available experimental data. Thus, the coupling function / cx a 2 certainly belongs to a 
class of models which lead to interesting scenarios. As a consequence, it seems natural 
to consider generalizations of the simple choice / oc a 2 . In fact, it will be proven in the 
following that the ansatz ( 1321) leads to a power-law for the spectrum of the magnetic 
field, the tilt n B being determined by the value of the parameter a (and n B = for 
a = 2) and, as a consequence, allows us to treat a simple but quite generic and general 
class of models. As mentioned before, the data are not yet very accurate and are still 
compatible with important variations of a around the preferred value a = 2. Given the 
present-day measurements, the corresponding range will be determined below. Then, 
given the ansatz fl32l) . we will reverse-engineer the spectrum and find the corresponding 
function /(</>) for a given model of inflation, taking into account the uncertainties on 
a. This will allow us to produce a complete class of models that satisfies the currently 
available experimental data. Then, in a second step, we will seek for particle physics 
realizations of these scenarios. 

Given Eq. fl32|) . the effective potential /"// can be expressed as 

^^>, (33) 

where we have defined 7 = a(l + f3). It is straightforward to integrate Eq. (121!) in terms 
of Bessel functions. The result reads 

A( V , k) = (kr]) 1/2 [d(A;) J 7 _ 1/2 (A^) + C 2 {k) J- 7+ i/ 2 (^)] , (34) 

where C\{k) and C*2(/c) are two scale-dependent coefficients which are fixed by the initial 
conditions. These ones are determined in the ultra-violet limit. Indeed, in the short- 
wavelength limit k/(aH) = —kn — > 00 , the vacuum reduces to the one in Minkowski 
space-time 

A(k, r,) - 4^ . (35) 



Notice that one can easily check that the normalization is correct. Indeed, the above 
relation means A(r], k) — > e~ tkv / (afV~2k) and taking into account the fact that A = A' /a, 
the previous mode function correctly reproduces the Wronskian given by Eq. ( 123]) (in 
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addition, [A] = [A] = M l / 2 as announced previously). As a consequence, the two 
arbitrary coefficients C\ and C2 are given by 

d(*0 = -C 2 (fc)e^- 1/2) , C 2 (k) = J — — ^ . (36) 

V 4/e COS(7T7j 

On large scales (compared to the Hubble radius), k/(aH) = —kr\ — > 0, using the 
asymptotic behavior of the Bessel functions, one obtains 

^ e t7r(l- 7 )/2 

+ 2-7+3/2 r(- 7 + 3/2) cos [^(1-7)]^ 2 (^) 1_7 . ( 37 ) 
where the amplitude of the two modes has been written such that the symmetry 
7 — ► 1 — 7 is manifest. Let us define the dimensionless function J 7 (5) by the following 
expression 

^ = 2 25 + 1 r 2 (5 + l/2) cos 2 (tt5) ' (38) 

with 8 = 7 if 7 < 1/2 and 5=1 — 7 if 7 > 1/2. This function is represented in 
Fig. [TJ If one assumes that the background space-time is almost de Sitter, i.e. (3 ~ —2, 
as indicated by the fact that the CMB power spectrum is almost scale- invariant, then 
one has 7 ~ —a. This function blows up at a = —1/2 which marks the frontier 
between the two branches. Otherwise, T{a) is smooth everywhere. In particular, 
the divergences due to the cosine at the denominator are exactly cancelled out by the 
singularities of the Euler function. For instance, if one considers the branch a > —1/2, 
this is best seen with the help of the equation, see Ref. [15], r (1/2 + S) cos (tt5) = 
T (1/2 — a) cos (ira) = tc/T (1/2 + a). At a = n/2 with n G N, the cosine vanishes but 
r (1/2 — a) = T(l/2 — n/2) blows up since the Euler function is singular for negative 
integer values of its argument. As proved by the above formula, the net result is finite 
since r (1/2 + n/2) is regular. 

Finally, the Fourier energy density of the magnetic field (in other words, the energy 
density stored at a given scale k) can be expressed as 

d k 3 1 / k x 25 

where, again, <5 = 7if7< 1/2 and 5=1 — 7 if 7 > 1/2 



= a?™? U I • (39) 



2.2. The Coupling Function 

In the previous section, we have derived the power spectrum of the magnetic field under 
the assumption that / oc a a . The goal of this section is to obtain the explicit form of 
f{4>) for various single field models of inflation. 

Let us start with power-law inflation since this is a situation where an exact solution 
is available. In this case, the inflaton potential is given by 

V(<j>) = V exp [-V2Ti (0 - <M] , (40) 
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Figure 1. Amplitude of the magnetic power spectrum (solid line) at the end of 
inflation, as given by the function J- (a) defined in Eq. (|38p . in terms of the index 
a characterizing the shape of the gauge coupling. An almost scale-invariant power 
spectrum, (3 ~ —2, has been assumed such that 7 = a(l + f3) ~ —a. The divergence 
at a = —1/2 signals the transition between the two branches of the spectrum. The 
amplitude of the electric power spectrum at the end of inflation is also displayed (dotted 
line). As shown in the following, the amplitude is given by the function Q(a) defined 
in Eq. (fTOO]) . 



where e\ is the first slow-roll parameter and is constant. Here, is dimensionless and 
denotes the vacuum expectation value measured is units of Mq- Then, in conformal 
time, the explicit solution of the Einstein equations reads 



|l+/9 



0(77) = o + y2^(l + /3)ln|77| 



(41) 



0(77) = £ c 

with (3 < — 2 and t\ = (2 + /3)/(l + (3). Therefore, we see that, in order to satisfy our 
requirement / oc a a , we must choose the coupling function such that 

a 



f((j>) oc exp 



(42) 



This is a priori interesting since this is precisely the shape postulated by Ratra in 
Ref . [30] . We will examine in Sec. H] whether we can design a well-motivated model 
which naturally gives this case. 

Let us now consider the case of a general potential. In this case, no exact solution 
is available but the slow-roll approximation can be used. Indeed, combining the two 
slow-roll equations of motion 



H 2 



-v{4>) 



dt 



1 d^ 



(43) 
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The integration of this equation is possible by quadrature. This provides us with the 
function a(<p) from which we trivially determine the form of the coupling function. One 
obtains 

' V((p) 



f{4>) oc exp 



—a 



V'{<p) 



dip 



(45) 



The above expression gives the form of the coupling function for any slow-roll model in 
terms of a single quadrature. In the following, we apply this general formula to various 
slow-roll scenarios. Let us start with the case of large field models, namely 

V{<f>) = M 4 p , (46) 

for which we find the following coupling function 

m ^(-gA . (47) 

Then, we can also consider the case of small-field and hybrid inflation where the potential 
can be written as 



V{$) = M 4 (l ± A0P) 



(4? 



where the (upper) plus sign corresponds to the hybrid case (p = 2 being the theoretically 
favored value) while the (lower) minus sign gives the small field case, and A is a constant. 
The corresponding coupling function reads 



f(4>) oc exp 



h2-p 



T- 



Xp(2 - p) 2p 
The case p = 2 must be treated separately. One obtains 

a0 2s 



exp 



(49) 



(50) 



The case of small field is particularly interesting. Indeed, in this case, <C m P1 , and the 
exponential term can be neglected. As a result, one finds 

/(0) oc 4> a/m ■ (51) 

The question is now to see whether one can find high-energy particle physics models 
which reproduce the previous phenomenological approach. This will be studied in detail 
in Sec. HJ 
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2. 3. The Magnetic Field during preheating 



In this subsection, we study the behavior of the magnetic field during the preheating 
stage. Previously, we have determined the shape of the magnetic power spectrum and, in 
the following, we will evolve it until present time using a simple adiabatic law. However, 
one must first check that preheating is not going to affect either the amplitude or 
even the spectral shape of dp B /dk. During preheating, at least if the potential can be 
approximated by a parabola, the inflaton field behaves according to 



cf>(t) 



\/3mt 



sin (mt) 



(52) 



In the following we work in terms of cosmic time since this is more convenient for our 
purpose. Then, the Fourier amplitude A(t, k) obeys the equation 



A(t,k)+ [3H + 2^)A(t,k)+ [% + H 2 + 2H f - + - 



f 



f 



A(t, k) = 0.(53) 



Let us now consider the new rescaled variable k) defined by A(t, k) = 

are led to introduce a new variable, slightly different from 
the one used before A(i], k), because we now work in terms of cosmic time rather than 
in conformal time. Using the above equation, it is straightforward to show that 

[ - 2 - { - H f - + \e 2 - ] W. />') = . (5 1) 



In order to go further, one needs to know the time-behavior of f(t) during preheating 
which requires the knowledge of the shape of the inflaton potential during this stage 
of evolution. In the simple approach considered here, this is the case only for chaotic 
inflation. Indeed, the shape of V((f) used in the power-law case (HO"]) or in the small 
field/hybrid case (148p is only valid during the slow-roll phase. Therefore, let us now 
concentrate on the case of chaotic inflation. In this case, using Eqs. (J4"T|) and (152]) the 
coupling function is given by 

4a . 9 

■ {mt) 



f(t) oc exp 



sin 



(55) 



3p(mt) 2 

From this expression, one can evaluate the time-dependence of the effective frequency 
in Eq. (I54p . For the term ///, straightforward calculations lead to 



/ 



4a 



+ 



3pt 2 

16a 2 
9p 2 m 2 t 4 
8a 



— - sin 2 (mt) + — sin(2mi) 



(mt) 



mt 



2 cos(mt) 



4 4 

— — sin 4 (mt) sin 3 (mt) + sin 2 (2mt) 

(mt) 2 mt 



cos(2mt) + O 



3pt 2 \mt 3 
As a consequence, Eq. (1541 can be approximated as 



— + -H 2 + 



8a 
Zpt 2 



cos(2mt) 



K(t, k) = 0, 



(56) 
(57) 

(58) 
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where we have retained the most and the least rapidly changing terms in the bracket in 
order to compare with the Mathieu equation, d 2 K/dz 2 + (a M — 2q M cos2,2) = 0. Then, 
defining z = mt — tt/2, Eq. ( 1581) is rewritten as 



d 2 N 

di 2 



+ 



k 2 H 2 8a 

+ n — 7T^ 777\~o COS 2z 



So we find 



a?m 2 Am 2 3p(z + 7r/2) 5 



k 2 H 2 4a 



K = . (59) 



~ M a 2 m 2 + 4m 2 ' ?M 3p(^ + vr/2) 2 ' (t>0) 

for a time scale much shorter than the Hubble time. Since we are interested in the long 
wave modes and H <C m, we find a M < 1. Then, the Mathieu equation has instability 
only for q M > 1. Since both a and p are constants of order unity with their typical 
values a = p = 2, we can conclude that the instability continues for less than one period 
of the field oscillation. The parametric resonance is therefore entirely negligible and we 
cannot expect any enhancement of the magnetic field during preheating, at least with 
the coupling function considered in this article. 

2.4- The Magnetic Field after Inflation 

After inflation and preheating, the Universe is full of charged particles and, therefore, 
the conductivity a c jumps up to a value much larger than the Hubble parameter almost 
abruptly. In this case, the model is described by the following action 

S [A,] = - J d'x^-g (^g^F^p - j%j , (61) 

where, now, the function f — 1 since the inflaton field has decayed and is no longer 
present. The current j M is defined by the formula 



jfi = Pe?V + CTcEf, . (62) 

In this expression, p e is the measurable charge density and a c is the scalar conductivity 
of the medium. The equation of motion now reads 

' d^{ y /=gg^^F uP )=f i (63) 



and the new stress-energy tensor can be expressed as 

T, v = -g^F^Ffr - ^"VW + W'/ - 2^ . (64) 
Then, the wave equation for the gauge field reads 

A(t, x) + (~ + a)j A { (t, x) - d j djMt, x) = , (65) 

In the large scale limit (hence neglecting the spatial derivatives of the gauge field) and 
in the large conductivity limit, cr c ^> H, the solution to the above equation reads 

Mt, x) = ^lM e -ct + D2i{x) ; (66) 

0"r 
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The exponential term will die away very quickly (with a characteristic time r = cr^ 1 ^> 
1/H) and, therefore, one obtains Aj(t, x) ~ D 2i (x). This implies that — and -B, is 
a constant in time. Let us now evaluate the energy density associated with the previous 
configuration. It is easy to see that the extra term coming from the current is given by 
p = a c A l Ai and therefore vanishes. At the end, only the magnetic part remains and 
reads 

P B = \^5 k \d 3 D 2 , - d e D 2j )(diD 2k - d k D 2t ) oc 1 . (67) 

Therefore, after inflation and until present times, one expects p B to scale as 1/a 4 
(independently of the era considered, ie radiation or matter J||. As a consequence, one 
deduces that the present magnetic energy density at a given scale L can be expressed 
as 

'■<* = °- £ > = t(* = *--* = T)Tfe) 4 - < 69 » 

Using the expression of p B at the end of inflation, one obtains, with the definition 
^ B ( k ) = P B ( Z = 0, A;)/p C ri, 

dfi fl 32 M/MVfiLfm (70) 



dlnfc 9 \ m Pl/ VPend/ \ a cnd/ V a 0-^l 

To go further, one must therefore evaluate the ratio a /a end . Clearly, this ratio depends 
on all the history of the Universe, in particular, on the process on reheating. It has been 
shown in Refs. [36J H7] that it can be expressed as 

^ L = U!i(»Lr"(3H!r 1/ \ (7i) 

Wend -ft 

where H ~ 2.1h x 10- 42 GeV, Q° rad h 2 ~ 4.3 x 10^ 5 , see Ref. p], and the parameter i?, 
which describes the reheating phase, is given by [l6l H7] 



1 - 3w re h 

In it — 



3(1 + w reh ) 



ln(^X h ) + 2( \t 3 3 ^ re \ ln ^ 



reh / 



(72) 



1 ln / 3 + 3 Wreh \ ^ 30 1 / 4 
1 - 3w re h V 5_ 3w reh / \/7r 
This parameter depends, a priori, on three quantities: the reheating temperature, T reri , 
the equation of state during reheating, w re h, and the energy density at the end of slow- 
roll inflation, p en d- However, if one considers a particular model then the number of 

| On large scales, one could worry that the Ohm's law is not applicable and, hence, the previous result 
not valid. However, even in this case, the conclusion that p B scales as 1/a 4 would be unchanged. 
Indeed, if the magnetic field is described by Eq. ([3]) but now with f = 1 since the inflaton field has 
decayed, then A = e~ tkri /y/2k and 

I I f + °° 

PsM = ^- 4 J q dfcfc3. (68) 

This integral diverges as M 4 , where M c is a cut-off as is usual for a vacuum contribution. But, as 
mentioned before, the important point is that it still scales as 1/a 4 . 



Generation of Large-Scale Magnetic Fields in Single-Field Inflation 



15 



free parameters can be reduced. For instance, in the case of large field models where 
V{4>) oc <f) p , the equation of state during reheating is known, namely 

Wreh = —T • (73) 



p + 2 

This leads to [3B| H7j 

3p 



In (# 1/4 T rch ) - lnp end 

v ' p — 4 



Moreover, the value of p en d can also be estimated. Indeed, the end of inflation is defined 
by the condition 

1 ^ 0e 2 nd /2 + ^(0e nd ) ' [ } 

which implies that 2 nd = V((p end ) or p cnd = W(<p cnA )/2. Since 4> e nd/m P1 = 1/(2^) 
for the potential V(^) = m 2 2 /2, one obtains p en d — 3m 2 mp 1 /(167r), the mass m being 
known, thanks to the COBE/WMAP normalization. Straightforward calculations lead 
to p cnd ~ 2160tt 2 /(2^ + l) 2 x Q 2 ns _ PS /T 2 ~ 1.03 x 10" 11 with Q ims - PS - 6x 10^ 6 , 
T ~ 2.7K and assuming that the number of e- folds between Hubble exit during inflation 
and the end of inflation is given by iV* = 50. We conclude that, in the case of a massive 
large field model, the parameter In R only depends on the reheating temperature. 

The situation can be more complicated for other potentials. For small field models, 
the shape of the potential during the slow-roll phase is different from the shape of the 
potential during reheating. As a consequence, the equation of state iy ren cannot be 
computed from the parameters describing V{4>) as it was the case for large field models. 
It remains a free parameter. In the case of power law inflation and of hybrid inflation, 
one needs a mechanism of instability in order to stop inflation. Then, even p en d remains 
a free parameter since it is linked to the value of the inflaton field where the waterfall 
behavior starts. Reheating proceeds in a direction perpendicular (in the field space) 
to the inflationary valley and, therefore, requires the full set (u> ren , T ren , p en d) to be 
phenomenologically described. 

Putting everything together one arrives at the following expression 



din A; 18?r 



dk) ' (76) 

/ i \ 25+4 

~ 2.4 xlO- 7 (2.28 xlO- 58 ) 8+2 F(5)R- 2S h 2S+2 ( . (77) 

\a H J 

This equation is one of the most important result of this article. It gives the amplitude, 
at a given scale, of the magnetic field today in terms of four parameters: 6, p em }, T Teh 

and W re h- 

In Fig. [2, we have represented the spectral index of the magnetic power spectrum, 
n B = 25 + 4 as a function of a as well as the quantity 5 (see the red dashed line). 
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a 



Figure 2. Spectral index n B (solid red line) of the magnetic power spectrum as a 
function of the index a, assuming, as before, a background expansion closed to de 
Sitter, i.e. (3 ~ —2. Scale invariance of the magnetic power spectrum corresponds to 
the values a — 2 and a = —3. The dashed red line represents the quantity d (see the 
text) as a function of a. The corresponding quantities for the electric field are also 
displayed, namely the spectral index n E (green dotted line) and the function t (green 
dotted-dashed line), see the definition of t in the text after Eq. (] 1 00[) and Eq. (jlOip . 
Scale-invariance of the magnetic power spectrum corresponds to a spectral index of 
the electric power spectrum given by n E = —2 for a — — 3 and n E = 2 for a = 2. Scale 
invariance of the electric power spectrum is realized for a — —2 or a = 3. 

The two branches of the spectrum are easily identified and meet at a = —1/2. Scale 
invariance is obtained for a = 2 or a = —3. For values of a such that —3 < a < 2, we 
have a positive index, that is to say a blue spectrum while for a < —3 and a > 2, we 
have a negative index i.e. a red spectrum. 

In the next section, we find the constraints put on a by the currently available data 
on the magnetic power spectrum. 

3. Constraining the coupling 

3.1. Limits on Cosmological Magnetic Fields 

Before comparing the above results with observation we summarize upper limits on 
cosmological magnetic fields from the following sources (see more detailed explanations 
in Refs. (HJHHJ). The first type of constraints comes from CMB anisotropy measurements. 
Indeed, homogeneous magnetic fields during the time of decoupling whose scales 
are larger than the horizon at that time cause the Universe to expand at different 
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rates in different directions. Since anisotropic expansion of this type distorts CMB, 
measurements of CMB angular power spectrum impose limits on the cosmological 
magnetic fields. Barrow, Ferreira, and Silk [49] carried out a statistical analysis based 
on the 4-year COBE data for angular anisotropy and derived the following limit on the 
primordial magnetic fields that are coherent on scale larger than the present horizon. 

^oLc < 5 x 1(T 9 G . (78) 



Moreover, in Ref. [50], a limit on primordial small-scale magnetic fields was obtained, 
also from CMB distortions. The constraint reads 

^oLc < 3 x 10" 8 G , (79) 

for coherence lengths from ~ 400pc and ~ 0.6Mpc. 

Another type of constraint comes from Big Bang Nucleosynthesis (BBN) since 
magnetic fields that existed during the BBN epoch would affect the expansion rate, 
reaction rates, and electron density. Taking all these effects into account in calculation of 
the element abundances, and then comparing the results with the observed abundances, 
one can set limits on the strength of the magnetic fields. The limits on homogeneous 
magnetic fields on the BBN horizon size ~ 1.4 x 10~ 4 Mpc are such that 

tf£Lc < 10" 6 G , (80) 

in terms of today's values [511 [52] . 

Rotation Measure (RM) observations also provide interesting limits. RM data 
for high-redshift sources can be used to constrain the large-scale magnetic fields. For 
example, Vallee [53] tested for an RM dipole in a sample of 309 galaxies and quasars. 
The galaxies in this sample extended to z ~ 3.6 though most of the objects were at 
z < 2. Vallee derived an upper limit of 6 x 10~ 10 (7?, e0 /10~ 7 cm~ 3 ) 1 G, where n e0 is the 
present mean density of thermal electrons, on the strength of uniform component of a 
cosmological magnetic field (n e0 < n b0 = (2.7±0.1) x 10~ 7 cm~ 3 , where n^o is the current 
density of baryons [54]). Let us also notice that Ref. [55] also quotes limits obtained 
from Faraday measurements, namely 

s£Lc < 6 x 10 " 9G . ( 81 ) 

for scales ~ 50Mpc and 

££Lc < 10" 8 G , (82) 

for coherence lengths corresponding to ~ G. 

Finally, one can also take into account a constraint which is of different nature. It 
is known [3] that the magnetic field on scales ~ IMpc must be 

a£Lc > 1(r22 ° > ( 83 ) 

in order to ignite efficiently the dynamo mechanism at the galactic scales, as was also 
mentioned in Sec. [TJ This constraint is clearly different from the above ones since it 
comes from a theoretical prejudice. However, it is interesting since this leads to a lower 
limit rather than upper bounds as discussed before. 
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Figure 3. Constraints on the amplitude of the present-day magnetic field at different 
scales. The bounds come from CMB measurements, Big Bang nucleosynthesis (BBN) 
and Faraday rotation (RM). The dynamo theoretical constraint is also plotted. The 
two solid blue lines represent the reddest spectra compatible with the CMB constraint 
at the Hubble scale in the case of a general model of inflation and in the case of 
large field models, respectively, while the dotted green line is the reddest spectrum 
compatible with the CMB constraint and the dynamo limit (valid in the case of a 
general model). Finally, the two dashed red lines are the bluest spectra (for a general 
model and for large field, respectively) compatible with the dynamo constraint. 



These limits can be translated directly into limits on Q B . Indeed, the energy density 
associated with the magnetic field reads p B = B 2 /2 (with our normalization of the gauge 
field kinetic term, see also Ref. [IE]) and this implies df2 s /dlnfc = 47r.B 2 /(3if ( 2 mp 1 ). 
Since lT 2 /2 ~ 1.9 x l(T 32 GeV 4 and IT = 10 4 G, one arrives at 



All the constraints on the present amplitude of the magnetic field are summarized in 



3.2. Constraining the Primordial Spectrum 

In this section, our main goal is to constrain the parameter 7 using the experimental 
data reviewed in the previous subsection. However, as already mentioned, the amplitude 
of the theoretical spectrum given by Eq. (176|) also depends on three other parameters, 
namely w re h, ^reh, Pend, that are not known precisely. Therefore, we first need to discuss 
the constraints existing on the values of those parameters. 





Fig. El 
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Let us start with the limits on the parameters R. First of all, there are limits 
on p end since, in order to preserve the success of the Big Bang Nucleosynthesis, one 
needs p en d > p nuc — 10~ 85 mp r On the other hand, one knows that the energy scale of 
inflation is constrained from above, namely H- lx ^/m vl < 1.3 x 10~ 5 , see Ref. jl6]. As a 
consequence, one must require that p enc ; < 10~ 10 m pl . Summarizing the constraints, one 
has 

- 187 < lnp end < -20. (85) 

Let us now consider the equation of state during reheating, u> re h- Reheating is, by 
definition, a phase of evolution where the scale factor does not accelerate which means 
that the strong and dominant energy conditions should be satisfied, that is to say 

- 2 < w rch < 1 . (86) 

Finally, one also gets constraints on p rch = g*TT 2 T^ eh /30 from the obvious requirement 
that reheating must proceed after the end of inflation and before the BBN, that is to 

say 

Pnuc < Prch < Pcnd • (87) 

Taking everything into account, one arrives at the following model- independent range 
of variation 

- In p nuc < In R < - — In p nuc + ^ In p end , (88) 
up to negligible factors depending on w re h only. This leads to 

- 47 < Initio. (89) 

In Ref. [46J, CMB data have been used to constrain the values of lni? for large, small, 
hybrid and running mass inflationary scenarios. Some very weak limits have been 
obtained in the case of small field models but, essentially, the currently available data 
are not yet accurate enough and the marginalized probabilities obtained are just cut 
at the edges of their priors. Therefore, in absence of any other constraints, one must 
consider the interval given by Eq. ( |89l) . This is what is assumed here for power-law 
inflation even if, strictly speaking, the analysis on In R should be redone for this model 
which was not explicitly considered in Ref. [16] . However, it seems extremely likely that 
no strong constraints would be obtained with the present-day data. In addition, let us 
recall that, in order to stop inflation, one must implement a mechanism of tachyonic 
instability where p en d is now viewed as a free parameter of the model. In this context, 
the constraint of Eq. (189]) seems to be particularly relevant. 

As already mentioned, the situation is slightly different for the massive large field 
model, p = 2. Indeed, in this case, the equation of state during reheating is known 
w re h = 0, as well as p cm } (at least approximatively, see above). In this case, one obtains 
the following constraint —20 < \nR < —6. In Fig. H] one has reproduced this range of 
variation together with the constraints on In R obtained by scanning general models of 
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Figure 4. Constraints obtained from the WMAP3 data on the parameter lni? at 
la- and 2a, see Ref. [IB]. The solid lines represent the variation of In R as a function 
of the energy density at the end of inflation for various reheating temperatures (or, 
equivalently, for various energy densities at the end of reheating). No observational 
constraint on T ro h is obtained in the case of large field models. 



the form, V{4>) oc P , for more details see Ref. |46j. This plot confirms that the reheating 
temperature remains a free quantity with the currently available CMB data. 

Let us now turn to the constraints on a. In the case where the spectral index n B is 
negative (red spectrum), the limiting constraint is the CMB constraint at scales of the 
order of the horizon today, see Fig. [3] So, one should require 



-10.2. (90) 



According to Fig. El n B < means either a > 2 or a < —3. Looking at the dashed red 
line, one sees that this implies 5 < in these two regimes and for this reason, one writes 
5 = —\S\. Then, using Eq. (1751) and working out Eq. ( |90l) . one obtains 

122.47 + log 10 (h*<m J din k\ CMB ) - log 10 T - 4 log 10 [k/ (a H )] 



\5\< 



57.92 + 0.86 In R - 2 log 10 [k/ {a H )} 

(91) 



where we have used h ~ 0.72. In fact, this inequality is not as simple as it may seem since 
T is also a function of 5. Therefore, given the functional form of J 7 , only a numerical 
calculation could allows us to derive the exact bound on 5. However, the impact of 
the term log 10 (jF) is limited and can be neglected in a first approximation. Moreover, 
since we seek an upper limit, it is clear from the above expression that the smallest 
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value of In R should be considered in order to reduce the denominator. According to the 
considerations presented before, this means lni? ~ —47 in general and InR ~ —26 for 
large field models. Then, using k/(aoH Q ) = 1 and log 10 (/i 2 df2 s /dln/c| CMB ) ~ —10.21, 
one arrives at \S\ < 6.1 and \S\ < 2.72 for large fields. The corresponding spectra are 
represented in Fig. [3J see the two solid blue lines. If a > 2, then 5 = —a and the two 
bounds respectively correspond to a < 6.1 and a < 2.72. On the other hand, if a < —3, 
then 5 = 1 + a. As a consequence, one obtains a > —7.1 and a > —3.72. However, one 
notices in Fig. [3] that the spectra with a = 6.1 and a = —7.1 violate the dynamo limit 
(this is not the case for large field models). Therefore, one should repeat the calculation 
and require the dynamo limit to be also satisfied. Straightforward calculations indicate 
that this amounts to a < 4.55 and a > —5.55 (with, this time, InR ~ —39.2). The 
corresponding spectrum is represented in Fig. [3] by the dotted green line. 

Let us now consider the case where we have a blue spectrum, n B > 0. According 
to Fig. [31 this means — 3 < a < 2. This time the situation could be more complicated 
because this range of values can correspond to positive or negative values of 5, see the 
dashed red line in Fig. [21 However, in practice, one remains in the regime where 5 < 
as will be checked below and, hence, 6 = —\S\. For blue spectra, the relevant constraint 
is the dynamo constraint which reads 
78.23 - log r 

1 1 50.68 + 0.86 In R ' v ' 

where we have used 

log 10 I K 2 cL ~ 



din A; 



37.61, (93) 

Dynamo / 

and log 10 [/c/(a -f^o)] ~ 3.62, the scale at which the dynamo constraint applies. Since 
we now seek a lower bound on one should use the largest value of InR, that is to 
say InR ~ 10 in general and InR ~ —6 for large field models. This respectively gives 
\5\ > 1.19 and \5\ > 1.55. This implies a > 1.19 and a > 1.55 in the case where 
a > —1/2 (and, hence, one verifies that, for these values of a, 5 is indeed negative, see 
Fig. [2]). In the other case, this means 5 = 1 + a < and, therefore, one has a < —2.19 
and a < —2.55. The corresponding spectra are represented in Fig. [3] by the two dashed 
red lines, the one with the smallest slope being the spectra in the case of large field 
models. 

Summarizing all the results obtained above, one obtains the following general 
constraints 

-5.55 < a< -2.19, 1.19 < a < 4.55 , (94) 

while, for the case of large field models, one arrives at 

-3.72 < a< -2.55, 1.55 < a < 2.73 . (95) 

This means that, if a is in the above ranges of values, then there exists at least one value 
of In R such that the corresponding spectrum is compatible with the currently available 
data. 
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3.3. Consistency and the Back-reaction Problem 



We now evaluate the electric field produced in the model under investigation. Our goal in 
this subsection is to check that the total amount of magnetic and electric energy density 
produced during inflation is not larger than the background energy density ~ Hf ni rri^ v 
Otherwise, clearly, the framework used here would not be consistent and would suffer 
from a serious back-reaction problem. Using the expression of the stress-energy tensor, 
see Eq. (12^1) . one obtains for the electric time-time component 

l/_ 2 
2 a 2 

from which straightforward calculations lead to the following expression for the vacuum 
energy density 



T 



EQ 



--^doAdoAj, 



(96) 



(T E %) 



1 f 



dkk' 



2tt 2 a 4 

or, in terms of electric energy density stored at the scale k 

_d_ f_)f_ \A(k,riY 
dk Pl ' 



2tt 2 a 4 



./ 



(97) 



(98) 



This equation should be compared to Eq. (|29|) . The solution for the Fourier amplitude 
of the vector potential has already been obtained previously in terms of Bessel functions. 
Using Eq. fl3ll) on large scales compared to the Hubble radius, that is to say k/(aH) = 
—kr] — > 0, one obtains 

. Prr p in~//2-in 



A' 



r 



I 



27+3/2 p( 7 + 3/2) C os(7T7) 



-i7r(3+ 7 )/2 



+ I,!, —, *-r\ F ; ^Hkv)" 1 , (99) 

2-7+1/2 r(_ 7 + 1/2) C os [-tt(1 + 7)] v U , V ) 

where, as was already done for the magnetic field, the amplitude of the two modes has 
been written such that the symmetry 7 — > —1 — 7 is manifest. Clearly, this symmetry 
is not exactly similar to the one obtained in the magnetic case. Let us now define the 
dimensionless function Q(l) by the following expression 



Q{f) = 



TT 



(100) 



2 2l+3 r 2 (i + 3/2) cos 2 (tu) 

with 1 = 7 + 1 if 7 < — 1/2 and 1 = —7 if 7 > —1/2. This definition is very similar 
to the definition (!38j) . The functions t(a) and G(a) are represented in Figs. [2] and [T] 
respectively. With the help of Q{l), we can write the Fourier energy density of the 
electric field. One obtains 



dk P ^ k) 



kr 



1 



k 



(101) 



27T 2 ' ' ' a 4 

This expression should be compared to the formula expressing the magnetic energy 
density at a given scale k, see Eq. ( 1391) . 
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We are now in a position where one can estimate when a back-reaction problem 
occurs. Clearly, the model is free of this difficulty when the following condition is 
satisfied 



din A; 



inf d In k 



<Pmt, (102) 
inf 



where, in the present context, the subscript "inf" means evaluated at the end of inflation. 
Indeed, as discussed previously, after the end of inflation, the conductivity jumps and, 
as a consequence, the gauge field becomes constant and, therefore, the electric field 
vanishes. Thus, if one checks that, at the end of inflation, the electric field can not 
cause a back-reaction problem, then we are guaranteed that the complete scenario is 
consistent. Let us also notice that the previous equation should be satisfied at any 
scales of astrophysical relevance today. Using the expression of the magnetic and electric 
energy densities derived before, the above relation amounts to 

0(0 



2tt 2 t 



I a H 



V a inf Hmi J \a H 

3 



o- \ 25+4 / , \ 2(5+4 

ao-no \ k 



<^-Hl f m 2 m . (103) 
ainfiiinf J \ ao-no / (N " 

In the following, for simplicity, we will not distinguish the energy density at the end 

of inflation from the energy density at which the scales of astrophysical interest today 

left the Hubble scale. Clearly, this is a good approximation during inflation (almost 

by definition). Under this assumption, using Eq. (ITTj) . one can estimate the ratio 

ao-ffo/(ainf-^inf) which enters the constraint ( 11031) . One obtains 

^^l.SlxHT 29 -. (104) 

Oinf-f^inf R 

Therefore, one sees that the constraint ( 11031) is quite difficult to analyze in full generality 
since, as already discussed, the quantity R depends on three parameters, namely the 
energy density at the end of inflation, the equation of state parameter during reheating 
and the reheating temperature, see also Eq. ( 1721) . In the following, instead of performing 
a systematic scanning of the parameter space, we just show that, in the vicinity of the 
"scale-invariant" values a = —3 and a = 2, there exist consistent models of inflation 
where the back-reaction problem does not occur. For this purpose, for simplicity, one 
first assumes instantaneous reheating, g*T* eh ~ 90iff nd m 2 >1 /(87r 3 ). Then, the expression 
of In R simplifies considerably and reads 

lni2~ ^lnp end . (105) 

Then it is straightforward to work out the constraint (11031) . For a = —3, using Eqs. (j 104ft 
and ( 11051) . one obtains 



Hinf \ , n _20 



2/3 



< Oil) x 10- 2U — — . (106) 
m P1 \a H J 

Since, in the present situation, the magnetic power spectrum is scale invariant, the above 
constraint comes almost entirely from the requirement that the electric energy density 
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Figure 5. Inflationary (solid line), magnetic (dashed line) and electric (dotted line) 
energy densities at the end of inflation in the case where a — —3, H m [/m m — 10 -22 
n B = and n E = —2. This situation corresponds to the limiting case treated in the 
text where there is no back-reaction problem, the electric energy density being always 
below the background energy density and where the dynamo constraint is still satisfied. 



should be less that the inflationary background energy density. Moreover, for a = —3, 
one has a red electric power spectrum with n E = —2, see Fig. [2j As a consequence, if 
the above constraint is satisfied at the Hubble scale, k/(aoH ) ~ 1, then it is satisfied 
at any scales of astrophysical relevance today. Therefore, for the case a ~ —3, one finds 
Hinf/m Pl < lCT 20 or p in f ~ lCT 41 m pl (let us recall that p nuc ~ lCr 85 m P1 ). 

However, we still have to check that the previous models satisfy the constraints 
of figure [31 Indeed, we have shown before that, for indices a in the ranges (|94"|) . there 
is always a value of lni?, compatible with the CMB data, such that the model is in 
agreement with the observations. However, we have not proven that the models are 
compatible for any values of lni? (this is clearly not the case) and, hence, since InR 
is now fixed (for instance with w rch = and H ini ~ 10 _20 m P1 , one has lni? ~ —21.77) 
one must verify that there is no problem with the data. In fact, it turns out that these 
models do not satisfy the dynamo constraint because the magnetic field is too strongly 
damped after inflation. If we want to solve this problem, one has to consider a different 
type of reheating era. Clearly, the most favorable situation is when u> re h ~ 1 since this 
is the situation where B scales the more slowly given that —1/3 < w rc h < 1, namely 
B oc t~ 2//3 and when the reheating stage is prolonged as much as possible, that is to 
say down to a reheating temperature of the order of few MeV's. In this case avoiding 
the back-reaction problem requires H in { < 10~ 22 m P1 which leads to lni? ~ —16.74. The 
corresponding situation is illustrated in Fig. [51 This time, for the previous value of In R, 
one can check that the branch a = —3 of the spectrum is still compatible with the 
dynamo constraint. Therefore, at least for this type of reheating period, this branch is 
still a viable alternative. 
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The case a ~ 2 can be worked out in the same manner. Straightforward calculations 
lead to 

^ < 0(1) x 10 58 ( — ^ 2 . (107) 
m P1 \a H J 

This time the electric spectral index is positive, n E = 2 and, therefore, the power 
spectrum is blue. As consequence, one should check that the consistency relation (11031) 
is satisfied at small scales, say k/(aoHo) ~ 10 10 . As is obvious from the above relation, 
this is very easily the case. 

The fact that back-reaction is always satisfied when a ~ 2 while it requires low 
scale inflation in the case a ~ — 3 can be easily understood, see also Ref. [39]. Indeed, 
from Eq. (|20|) . one has Ai(r],x) ~ 6i\A(r], k) ~ A(r],k)/ f(r]), where we have used 
the fact that the covariant polarization vector contains a factor a(rj), see Eq. ( TT8l) . 
Then, from Eq. ff2TT) . one has the two following modes on larges scales A oc f(rj) and 
A oc fiji) p dr// 2 (r). Clearly, from the above equations, the first mode leads to 
Ai ~ Const, and, therefore, to a vanishing electric field, while the second remains time- 
dependent and can be responsible for the production of a strong electric field. Indeed, 
it is easy to see that, in this last case, Ai(r], x) oc a -1 ~ 2a , that is to say Ai(r],x) oc a 5 
for a ~ —3 and Ai(rj,x) oc a~ 5 for a ~ 2. Therefore, this mode is sub-dominant for 
a ~ 2 but is dominant for the other branch of the spectrum a ~ —3. In this case, its 
strong time-dependence oc a 5 produces a strong electric field and this explains why the 
back-reaction problem can be severe in this situation. This problem can be avoided if 
the production of the electric field is limited as it happens to be the case if inflation 
proceeds at low scale. 



4. Particle Physics Models 

4-1. Power Law Inflation 

In this section, we finally turn to the model-building issue. More precisely, the question 
addressed in this subsection is whether there is a natural high energy model with an 
exponential inflaton potential and an exponential gauge coupling, as was obtained in 
Sec. 12.21 We will consider the other models in the next subsections. In order to study 
this issue, let us consider that the inflaton field is a modulus field T in string theory. 
Then, its Kahler and super potential are given by 

M 3 

K = -pin (T + T f ) , W = T n , (108) 

v ' n 

the natural choice being p = 3 (no-scale). Moreover, in string theory, a natural choice 

is f(T) = T q with the preferred choice q = 1 [H] . Therefore, at first sight, this does 

not give the required / function. However, one has to remember that the field T is not 

canonically normalized because the Kahler potential is non-trivial. The canonical field 

<p is given by 

T = e -v /2 ^. (109) 
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As a consequence, one obtains a coupling function with the required exponential shape, 
namely 

f(4>) oc e - q V^* . (110) 

Moreover, the F-term potential is also easily computed and one also obtains the required 
shape, namely an exponential dependence. The explicit expression reads 



V(cf>) 



M 6 (n - 3) 



Qn 



exp 



(2n-3)W-0 



(HI) 



If we compare the two previous formulas to Eqs. (142]) and (140]) . then one obtains the 
correct power-law model for the choice 

"-KWD- (112) 

In particular, we see that the branch a ~ 2 leads to a negative power for q while the 
branch a ~ —3 gives q = 3v / 3/(2y / ej") ^> 1. This is already an indication that it is 
difficult to reconcile the branch a ~ 2 with a sensible model building. In fact, in the 
present context, even the a ~ —3 branch cannot lead to the string inspired favored 
model with q — 1. But, the lethal flaw of the above model is that the required value of 
n is less than 3 leading to a negative inflaton potential. 

The next question is whether the mechanism of assisted inflation [56] could solve 
this situation. Let us envisage a situation where we have N moduli such that 

N N 



i=l 



W 



- E — T "* 



i=i 



rii 



(113) 



the no-scale structure being preserved if Y^i=iPi = 3- As before, the field T, are not 



canonically normalized, the normalized field 



being given by 



(114) 



Straightforward manipulations lead to the following potential 



V 



M 6 



exp 



N \ [ N N 1 

£V2^J -££- 

J=l / L j=l k=l K 




exp — rij 



N 




(115) 



/>/,W —<.>i,- I + y~] — exp 

If one considers the case of one field (which implies p — 3), then one can check that 
Eq. f illip is correctly recovered. Then, let us consider that we have N fields such that 
Pi = 3/N and rij = n for any index i. Then, the potential becomes 

M 6 N 2 (n -3) 



This implies that 



n 



3N 



6n 



AT3/2 



exp 



(2n - 3N) 



2N 



(116) 



(117) 
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We see that N > 2 allows us to have a positive potential and, therefore, represents 
a meaningful model. Of course, the fact that the power of the superpotential is not 
an integer is not a nice feature and, for this reason, the model remains unsatisfactory. 
What about the gauge coupling? In presence of several superfields, the form of / that 
should be chosen remains unclear. If we postulate that it depends on one superfield 
only (say the field then one gets 



In this case, and contrary to the single field case, q is not necessarily very large (as 
one might think since it scales as the inverse of the slow-roll parameter) because pi can 
be a small number. In the case envisaged above, this is generically the case because 
Pi = 3/N (remember that only the condition YliLiPi = 3 must be fulfilled in order to 
preserve the no-scale structure). Moreover, one could also imagine a case where all the 
coefficients pi are not equal. If one of the pi is relatively small, then it could compensate 
the inverse of the slow-roll parameter in order to give a number of order one. For the 
branch a ~ 2, q is negative and we are not aware of any well-motivated SUGRA model 
where this happens. Therefore, the most interesting case seems to use the other branch 
of the spectrum with a ~ —3. Then, one obtains 



with | a | ~ 3 and the favored model q — 1 is obtained for 



Of course N should be an integer and this can be the case if a is not exactly —3. We 
have seen that this situation is perfectly compatible with the data, see the previous 
section. Let us give an example for the sake of illustration. We can take (3 ~ —2.025 
leading the spectral index n s = 0.95. Then, the value a = —3.00135 gives N = 277 
fields in the model. 

A last remark is in order here. We have seen that the branch a ~ —3 is a viable 
solution only if the scale of inflation is low and the reheating stage quite long. Here, we 
can certainly design a model where this is the case since, in order to stop inflation, one 
needs to rely on a mechanism of tachyonic instability as in the case of hybrid inflation. 
In this situation, the energy density at the end of inflation remains a free parameter and 
can be chosen such that it leads to the correct scale of inflation. Moreover, the details 
of the reheating will depend of the shape of the potential in the direction perpendicular 
to the inflationary valley. 

4-2. Large Field Inflation 




(118) 




(119) 



One can generate large-scale magnetic field in large field inflation models if gauge kinetic 
function has appropriate dependence on the inflaton field given by Eq. fj47|) . As the 
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model of inflation we adopt a chaotic inflation in supergravity based on a shift symmetry 
$ — > $ + iC where C is a real number. Specifically we adopt a model proposed in 
Ref. [57] with two chiral superfields $ and X. In this model by virtue of this symmetry 
and its soft breaking, the Kahler potential and superpotential are given by 

K = + + IX* + --- , W = mX<5>, (121) 

where m is a constant corresponding to the inflaton mass. The scalar Lagrangian 
possesses standard kinetic terms and is given by 

^scalar = -dpQW ~ d^Xd^X* - V($, X) , (122) 

where the potential reads 

V{$,X) =m 2 e x ||$| 2 (l + |X| 4 ) 

+ \X\ 2 [1 - |$| 2 + ($ + $*) 2 (1 + |$| 2 )] | . 



(123) 



Decomposing the complex scalar field $ into two real scalar fields as $ = (a + icj)) / y/2 
and using the fact that only <p can have a large initial value beyond the Planck scale 
without any exponential barrier due to the Kahler potential, we find that the Universe 
evolves as in the chaotic inflation model with a massive scalar field |58j, namely 

Anflaton = " V(<j>) , V{4>) = ^m 2 <f) 2 . (124) 

Therefore, as announced, one has obtained the potential (j4Ul) with p — 2. Here, we 
recall that all models such that p > 3.1 are now excluded at 95% confidence level by 
the WMAP3 data, see Ref. [16]. 

We now consider the vector part of the Lagrangian. Specifically we take the 
following Lagrangian for the U(l) gauge field 

£ EM = -^/»F F ^, (125) 

with 

/(*) ee / e Q * 2 , (126) 

where a is a constant. Note that the shift symmetry of the system is preserved if 
Ffj, v = 0. As argued above, since a rapidly relaxes to the origin at the onset of inflation, 
the gauge kinetic function then reduces to 

»/(*) = /oe"^ 2/2 = f\<P) , (127) 

and it reproduces the coupling function given by Eq. (|47p when p — 2. However, and 
contrary to the case of power-law inflation, the choice (11261) can no longer be justified 
from string theory, at least for the simple favored case mentioned before. 
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4-3. Small Field Inflation 

Generation of large-scale magnetic field is also possible in small field inflation models. 
Here we adopt a model proposed in Ref. [59] which does not require fine tuning of 
the initial conditions. This model is based on the following Kahler potential and the 
superpotential 

K = -($ + <F) 2 + XX* , W = v 2 X (1 - g<t> 2 ) . (128) 

Symmetry argument to obtain the above expression has been fully described in Ref. [59] 
and we do not repeat it here. In this model the real part of the scalar component 
$ = (0 + ix)/y/2 plays the role of the inflaton for the small field inflation while 
its imaginary part serves as that for chaotic inflation which occurs beforehand. The 
quantity X is another chiral superfield introduced to yield an appropriate scalar 
potential. The scalar part of the Lagrangian then reads 

with the potential V(<p,x,X) given by the following expression 



£(0, x , X) = -d^<j> + -dpx&x + dpXdPX* - V((j>, X , X) , (129) 



V(cf>,x,X)=v*eW + * 

x (1- \X\ 2 + \X\ A ) + \X\ 2 



i 9 i2\ 2 i 2 /i i 9 ,2 i 9 2 

l ~2 < t ) ) +9X [1 + +4X 



2< ? Y + 2( 5 -l)V 



+ 2g(g + 1) <p 2 x 2 + 2g(g- 1) 4 + ^ (0 2 + X 2 ) 2 



(.130) 



Since x is f ree from exponential rise as in the model discussed in the previous subsection 
it can naturally induce chaotic inflation. Meanwhile and X settle to the origin apart 
from quantum fluctuations. After chaotic inflation is still close to the origin and then 
induces new inflation. In this regime with x — \X\ = and \<p\ <C 1, the scalar potential 
reads 

Vtf) =v'(l- C -<P 2 ) , (131) 

with c = 2(g — 1). Thus, if g > 1, (j) rolls down slowly toward the vacuum expectation 
value 0o = y2~Jg and small field inflation takes place. The potential (I13ip has the same 
form as (1481) . 

Finally, using Eq. ( IBTl) . the desired form of the gauge kinetic function is realized if 
we take 

/>) = /o$ 2Q/c . (132) 

The branch a ~ 2 is the one to be used in order to obtain positive q. The favored string 
inspired model q = 1 is obtained for c ~ 4 or g ~ 3. However, unfortunately, this value 
is not compatible with the CMB constraints. Indeed, in the case of small field models, 
the first slow-roll parameter e\ is exponentially small while the second one is given by 
e 2 = 2c, see Ref. [46]. The CMB constraints on e 2 are -0.07 < e 2 < 0.07 at 95% CL, 
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see Ref. |46j . which implies —0.035 < c < 0.035. This means that q > 57, a value not 
compatible with the string inspired value and too large to be considered as natural. 
However, at the same time, the model makes use of the branch where the back-reaction 
problem is not severe at all. In this case, the characteristics of the required reheating 
stage seems easy to satisfy. 

5. Conclusion 

In this section, we recap our main results and discuss issues that should be studied 
in future works. In the present paper, we have studied the generation of large-scale 
magnetic fields in inflationary cosmology, breaking the conformal invariance of the 
electromagnetic field by introducing a coupling with the inflaton field as it is generic in 
a supergravity framework. We have determined the form of the coupling necessary in 
order to produce large-scale magnetic fields with the required strengths on the relevant 
astrophysical scales. We have shown that the scale-dependent magnetic energy density 
possess two branches. Among these two branches, one (a ~ —3) seems to lead to sensible 
model building in the case of power-law inflation but only at the expense of having a 
low energy inflation scale and a long reheating stage. The other branch (a ~ 2) is useful 
in the case of small field models and this case appears to be relevant since this does 
not require to fine-tune the reheating epoch. However, the model building condition 
q > 57 cannot be justified in a string inspired model and seems artificial. Finally, in the 
context of large field models, we have shown that the required coupling is quite difficult 
to justify. 

Determining exactly all the consistent models in the framework envisaged here is 
a non-trivial issue because the parameter space is large (at least four-dimensional), in 
particular when a general reheating stage is considered. Here, we have just proven that 
consistent models exist. However, it would certainly be interesting to systematically 
explore the parameter space in order to have a more accurate idea of whether the 
consistent models are just peculiar or, on the contrary, quite generic. Another interesting 
avenue for the future is clearly the model building issue. Here, we have noticed that 
the branch a ~ 2 can be used only in the context of small field models. It would be 
interesting to find other models where this can also been done. We hope that these 
issues will be addressed in the near future. 
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